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ABSTRACT 

Extended 21cm absorption regions (dubbed "21cm absorption halos") around first 
galaxies at z ~ 30 are likely the first distinctive structures accessible to radio obser- 
vations. Though the radio array capable of detecting and resolving them must have 
~ 200 km^ total collecting area, given the great impact of such detections to the un- 
derstanding of the reionization process and cosmology, such radio survey would be 
extremely profitable. As an example, we point out a potentially useful byproduct of 
such survey. The resolved 21cm absorption "halos", likely close to spherical, can serve 
as (almost) ideal sources for measuring the cosmic shear and mapping the matter 
distribution to z 30. We investigate the expected lensing signal and consider a va- 
riety of noise contributions on the shear measurement. We find that S/N 1 can 
be achieved for individual "halos". Given millions of 21cm absorption "halos" across 
the sky, the total S/N will be comparable to traditional shear measurement of ^10^ 
galaxies at z ~ 1. 
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1 INTRODUCTION 



Gravitational lensing (se e. 



revi ews of 

iBartelmann fc Schneider! I2OOII : iRefregieil |2003| ) directly 
probes the matter distribution of the universe and is 
becoming one of the most important probes of cosmology. 
Gravitational lensing distorts galaxy shape {cosmic shear), 
changes galaxy number density {cosmic magnification) 
and induces mode-coupling in cosmic backgrounds. These 
observable lensing effects enable several powerful methods 
to extract cosmological information. Cosm ic mic r owave 
background (CMB) l ensin g tSeliak fc Zaldarriagal Il999l: 



Zaldarriaga fc SeliakI 1 19991 : IHu 



2001 



Hu fc Okamotd 



20021), 21cm backgr o und le n sing dCo orav' '2004'; 'Pen' '2004|; 



Zahn fc Zaldarriagal I2OO6I : iMandel fc Zaldarriaga 200® 



and cos mic magnification measure ments of 21cm emitting 
galaxies (|Zhang fc Penll2005l .[ 20061 ) are expected to achieve 
high signal-to-noise ratio (S/N) in the near future. On the 
other hand, cosmic shea r measurements have achieved high 
S/N (e.g., [jj ■ 



vis et all [2OO6; Van WacrbokelitaD l2005l : 



iHoekstra et al. 20061) and will be improved significantly by 



several ongoing or upcoming large lensing surveys. 

Traditional cosmic shear measurements are fundamen- 
taUy limited by intrinsic ellipticities of galaxies. Intrinsic 



ellipticities have a dispersion ~ 30%, much larger than the 
typical 1% lensing signal for source galaxies at 2 = 1. Even 
in the best case that ellipticities of galaxies do not correlate, 
one still needs to average over several hundred galaxies to 
achieve S/N of one. Furthermore, source galaxies at zs > 3 
are difficult to detect opticaUy, so it is hard to map the mat- 
ter distribution at redshifts beyond 2 in optical lensing sur- 
veys. These two intrinsic problems can be over come by th e 
21cm absorption regions around first galaxies ('Cen"2006!). 
These regions are of arc-minute size, or ~ Mpc in comoving 
scale. Since their sizes are much larger than the nonlinear 
scale at corresponding redshifts, they are believed to be very 
close to spherical (Ccn 2006). In addition, since they lie at 
redshifts z ~ 25, the lensing signal is strong. These two in- 
trinsic advantages make these regions nearly ideal targets 
for shear measurement. Since these regions have distinc- 
tive structures, they are dubbed "21cm absorption halos" or 
"21cm halos". We caution the readers that these 21cm "ha- 
los" are not virialized dark matter halos. In fact, the density 
fluctuations in these "halos" are well in the linear regime. 
In this paper, we study the requirement and application of 
the shear measurement of these 21cm "halos". 

Observations of 21cm absorption "halos" are challeng- 
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ing and are way beyond the capability of the planned square 
kilometer array However, mission capable of resolv- 

ing these "halos" will be extremely profitable. The 21cm 
brightness temperature can be simultaneously measured to 
high accuracy over a large fraction of the sky and a wide 
redshift range, without extra cost. This will allow the mea- 
surement of related statistics, such as the brightness temper- 
ature power spectrum, bispectrum, etc., to unprecedented 
accuracy and significantly improve our understanding on 
reionization process and cosmology at high redshifts. The 
detection of 21cm absorption "h alos" will h ave potentially 
great impact on cosmology too (|Cenl [20061 ). Furthermore, 
the same mission can detect even smaller 21cm absorption 
regions around first stars, the so called Lyman-a spheres 
(because it is mainly the Lyman-a scattering that cou- 
ples the 21cm spin temperature to the kinetic tempera- 
ture, leading to the 21cm absorption against the CMB, 
IChen fc MiraldaJ ^scudc 200Q), through strong lensing of 
galaxy clusters I Li et al...2007i ). Given these potentially pow- 
erful applications, such mission deserve detailed investiga- 
tion and observation efforts. In this paper, we discuss a sur- 
prisingly valuable byproduct of such mission: mapping the 
matter distribution to z ~ 25 through weak lensing of 21cm 
absorption "halos". 

As we show in this paper, the ellipticity induced by var- 
ious contaminations is very likely to be controlled to below 
10% for 21cm absorption regions considered here and we 
expect that shear measured from each such region would 
achieve S/N~ 1. At z ~ 25 with Az = 2, for a lower mass 
cut ~ 7 X lO^MO, the total number of these 21cm absorp- 
tion "halos" is ~ 10*^ (Fig.[T)|. The total S/N will be equal 
to that with ~ 10^ z ~ 1 source galaxies and thus compara- 
ble to what would be provided by planned ambitious cosmic 
shear surveys such as the Large Synoptic Survey Telescopfl 

Moreover, the lensing applications of 21cm absorption 
"halos" will fill the gap of source redshifts between tradi- 
tional cosmic shear measurements (source redshifts zs ^ 3), 
cosmic magnification of 21cm emitting galaxies {zs ~ 1-6) 
and CMB lensing {zs ~ 1100). They provide independent 
checks for 21cm background lensing {zs ~ 10 — 30). 

Furthermore, lensing measurement to such high red- 
shifts will provide a direct way of CMB delensing and 
improve the me asurement of primordial in flationary grav- 
itational waves l|Sigurdson fc Cooravl[2005l ). Lensing recon- 
struction from 21cm absorption "halos" will significantly im- 
prove our understanding of matter distribution to z ~ 25 
and the nature of dark matter, dark energy, gravity and in- 
fiation. 

The paper is structured as follows. After a brief intro- 
duction of the 21cm absorption "halos" of first galaxies in 
§ 2, we discuss their detection prospect in § 3. Then, in § 4 we 
define the mean cosmic shear to be measured for these "ha- 
los". It would be ideal for cosmic shear measurement if the 
21cm "halos" are intrinsically spherical. Therefore, in § 5, 
by investigating several sources that may induce ellipticity 
of these "halos", we estimate the expected ellipticity and 
show that these regions are quite close to spherical. Given 
the expected noise, in § 6 we show the range of lensing power 
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spectrum that can be precisely measured. Finally, we sum- 
marize and discuss the main results. Throughout the paper, 
we assume a spatially flat ACDM cosmology with matter 
density parameter flm — 0.26, baryon density parameter 
Qb — 0.044, primordial power index ns = 0.95, Hubble con- 
stant HO = 100^ = 72km s~^ Mpc~^, and a fluctuation 
amplitude crS — 0.77 on 8/i~^Mpc scales, which is consistent 
with the Wilkinson Mic rowave A nisotropy Probe (WMAP) 
three- year results (S pergel et al.|[200Q ). For some plots, we 
also show cases with o"8 = 0.9. 



2 21CM ABSORPTION "HALOS" OF FIRST 
GALAXIES 

Since the standard cold dark matter (CD M) cosmologi- 
cal m odel is by now accurately determined (|Spergel et al.l 
l2006h . one could make testable predictions for the first gen- 
eration of galaxies, expected to form within halos of mass 
M = 10^— lO^M© at z = 20—40 (Fig.©. The first galax- 
ies are expected to be small and faint in terms of stellar 
optical light. The very first galaxies may have masses of 
~ 10^ M0 (Peebles 1980), limited by Jeans mass and molec- 
ular hydrogen cooling. However, star formation in minihalos 
with molecular cooling might be quenched, when hydrogen 
molecules are destroyed by Lyman- Werner photons. In any 
case, because of the expected low star formation efficiency 
in the minihalos, their 21 cm signals will be too weak to 
be relevan t even if they form stars. However, as shown in 
ICenI (|2006l ), with the coupling of the 21cm spin temperature 
to Lyman-a scattering, some of the large first galaxies with 
total mass M ^ lO'^-^M© may each display an extended 
hydrogen 21cm absorption "halo" against the cosmic mi- 
crowave background with a brightness temperature decre- 
ment of 5T = -(100—150) mK at a radius of 0.3—3.0 Mpc 
(comoving), corresponding to an angular size of 10 — 100 arc- 
seconds. These large galaxies can cool by atomic processes. 

The detection of these radio "halos" around first galax- 
ies will be extremely profitable, although admittedly diffi- 
cult. If and when a 21cm radio survey of first galaxies is 
carried out, some fundamental applications for cosmology 
and galaxy formation may be launched, as discussed in Cen 
(2006). Here, we recapitulate them. First, it may yield di- 
rect information on star formation physics in first galaxies. 
Second, it could provide a unique and sensitive probe of 
small-scale power in the standard cosmological model hence 
physics of dark matter and infiation, by being able to, for 
example, constrain the primordial power index ns to an ac- 
curacy of Ans = 0.01 at a high confidence level. Constraints 
on the nature of dark matter particles, i.e., mass or tem- 
perature, or running of index could be still tighter. Third, 
clustering of galaxies that may be computed with such a sur- 
vey will provide an independent set of characterizations of 
potentially interesting features on large scales in the power 
spectrum including the baryonic oscillations, which may be 
compared to local measurements l|Eisenstein et al.]|2005l ) to 
shed light on gravitational growth and other involved pro- 
cesses from z = 30 to z; = 0. Finally, the 21cm absorp- 
tion halos are expected to be highly spherical and trace the 
Hubble fiow faithfully, and thus are ide al systems for an ap- 
plica tion of the Alcock-Paczyriski test (|Alcock fc Paczvhskil 
I1979I ). Exceedingly accurate determinations of key cosmo- 
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Figure 1. Solid curves are cumulative halo mass functions at 
z =30, 2 5, and 20 (b ottom up), predicted by the Sheth-Tormen 
formula llSheth fc T ormen 200:^). Dotted curves are products of 
the star formation efficiency c* and the fraction /coll of mat- 
ter that has collapsed to halos where stars has formed at the 
three redshifts. The horizontal line represent a threshold value 
of c * /coll, below which the heating of the IGM by the X-ray 
background is not sign ificant, with a temperature increment less 
than ~2K l|Cenll200d ). This figure is a reproduction of Fig. 3 in 
ICenI 1I2OO6I) using the more accurate Sheth-Tormen mass function, 
assuming c* = 0.2 (0.001) for large halos (minihalos), and adopt- 
ing cos mological paramete rs consistent with WMAP three-year 
results jSpergel et al.ll2006l) . Om = 0.26, UA = 0.74, Qb = 0.044, 
h = 0.72, a8 = 0.77 and ns = 0.95. We focus on halos with mass 
greater than ~ 7x lO'^M©, which could harbor first galaxies capa- 
ble of prod ucing extended 21cm absorption regions against CMB 
llCenll2006D . 



logical parameters, in particular, the equation of state of 
the dark energy, may be finally realized. As an example, 
it does not seem excessively difficult to determine the dark 
energy equation of state w to an accuracy of Aw ~ 0.01, 
if flm has been determined to a high accuracy by different 
means. If achieved, it may have profound ramifications per- 
taining dark energy and fundamental particle physics (e.g., 
lUpadhve et al.l 120051 ) . This last property may also make 
them ideal sources for weak gravitational lensing measure- 
ments, as investigated here. 



3 DETECTABILITY 

The resolution and r.m.s. noise of a radio survey depend on 
the detailed array configuration. Without loss of generality, 
we focus on a configuration oi N x N arrays homogeneously 
distributed over a squ are with side length L. The resolu- 
tion is 9p = yj 7ra^/4A/L, where a ~ 1.2 for this specific 
configuration and A = 0.21(1 + z)xa is the redshifted wave- 
length of the 21cm line. Foreground contaminations at the 



low frequencies of concern here are overwhelming. The domi- 
nant one is galactic synchrotron emission, which scales as 
with /3 ~— 2.55. For concreteness, we adopt the brightness 
temperature at i/~54MHz (redshifted frequency of the 21cm 
lines at z = 25) as T^'^^ ' b ~ 3 000 K. This number is consis- 
tent w ith lkeshet et all (|2004l ) and lChen fc Miralda-Escudel 
(|2006l ). However, brightness temperature estimation at these 
low frequencies is qu ite uncertain (e.g., r°^"6 adopted by 
I Bowman et al.l I2OO6I is a factor of 2 higher, if scaled to 
z = 25.). A factor of 2 increase in T^^'^b would require a 
factor of 4 increase in integration time or a factor of 2 in- 
crease in the total collecting area, in order to achieve the 
same S/N. 

Since the foreground is highly smooth in frequency 
space and the 21cm signal has line features, showing as sharp 
fluctuations in the spectrum, the mean foreground contami- 
nation can be efficiently removed pixel by pixel ( Wang ct al] 
I2OO6I ). The residual noise per resolution pixel (with area 6p ) 
caused by photon number fluctuations in f oreground and in- 
strum ental noise can be worked out to be (jThompson et all 
I2OOII) 



aT 



Tsys 4:V2 
y/Aut na"^ f cover 

Tsys / Af 



(1) 



10 mK 



3000K 



\ 300kHz lyear J \ 



/cover \ " 



Here, /cover = Atotal/L^ = N^Adish/L^ and Adish is 
the collecting area of each dish. The bandwidth Ai/ is re- 
lated to the comoving separation r by Au ~ 46kHz [26/(1 -I- 
z)]^/^[r//i-^Mpc]. The typical bandwidth Av = 300kHz 
corresponds to ~ 6/i~^Mpc (comoving) at z ~ 25, which 
is about the diameter of the 21cm absorption "halos" (see 
Fig. [2)). The integration time t per pixel is also the integra- 
tion time per field of view (FOV). The system temperature 
Tsys ~ TV^'^ in our case. We note that the noise estimate 
is much larger than that from the Poisson noise of photons. 
The reason is simple — in the Rayleigh- Jeans part of the 
emission spectrum, the mean occupation number n of pho- 
tons per quantum state is much greater than unity and the 
variance, which is n(n -1-1), is super- Poisson. 

With assumptions on the star formartion properties in 
the central sources, the brightness temperature pr ofile of the 
21cm absorption "halos" can be computed (see ^Cen! !2 006l 
for details). At z = 25, 21cm absorbing regions generated 
by first galaxies in 1O*M0 halos have angular size of sev- 
eral arc-minutes with signal < — 50mK (top panels. Fig. 
To reliably measure the shape of these "halos" , a resolution 
of 9p ~ 0.3' and a system temperature noise aT ~ lOmK 
are required. For Av = 300kHz and /cover = 14%, this 
translates to L ~ 40 km, j4total ~ 200km'^, and an inte- 
gration time of one year per FOV. Detecting such regions 
at 5-a level without resolving them is relatively easy, which 
requires L ~ 4 km and ^total ~ 2 km^. The FOV at these 
wavelengths is typically ~ 10^ deg'^. To cover the whole sky, 
the total integration time of the order of half a century is 
required. Clearly, resolving these 21cm "halos" is not doable 
in the near future without major technological breakthrough 
or significant improvement in financial resources. Still, the- 
oretical merit of this study remains and could be useful for 
mid and far future proposals on dark age science. 

The signal and size of the 21cm absorption "halos" do 
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Figure 2. Top left panel: the brightness temperature of a 21cm 
absorption "halo" as a function of radius (in comoving unit), 
around a first galaxy in a 10^ Af© halo at 2: = 30, 25, 20, respec- 
tively. The star formation efficiency c* = 0.2 is adopted. Top right 
panel: the brightness temperature averaged over ±3/i~^Mpc (co- 
moving) along the line of sight. At these redshifts, 1 arcminute 
corresponds to ~ 2.4/i~^Mpc (comoving). Bottom left panel: el- 
lipticities caused by random system noise (with an r.m.s. aT per 
resolution pixel). See the text for details. Bottom right panel: the 
lensing signal, k or 7, smoothed over an area . The arrow de- 
notes the size of the 21cm absorption regions used for shear mea- 
surement. The dashed line is roughly the expected noise level. 



not have strong dependence on the redshift (Fig. [2} , until 
X-ray heating becomes important. However, towards higher 
redshifts (e.g., z — 30), thus lower frequency for observation, 
the foreground increases significantly. Detecting and resolv- 
ing these regions at z = 30 would require more ambitious 
surveys. Furthermore, the number of halos drops by a factor 
of ~ 100 from z = 25 to 2: = 30 (Fig.[T]), which substantially 
reduces the significance of shear measurement at z = 30. 
These factors make the measurement of shear for 2: = 30 
"halos" very unlikely. Towards lower redshifts, the shear 
measurement is also unlikely, but for a different reason. One 
necessary condition for the existence of the 21cm absorp- 
tion feature around the first galaxies is that the intergalac- 
tic medium (IGM) is not significantly h eated by th e X-ray 
background. We recast the arguments in [Ceiil (I2OO6I) by us- 
ing t he more accurate halo mass function (|Sheth fc TormenI 
|2G02| ) and adopting cosmol ogical parameters co nsistent with 
WMAP three-year results (|Spergel et al.ll2006l ). The resuhs 
are shown in Figure [TJ and we find that the above condition 
can be satisfied at z > 23. Therefore, we will focus our study 
for z = 25. 



4 COSMIC SHEAR 

To estimate the feasibility of shear measurement from 21cm 
absorption "halos" , we make several simplifications. We as- 
sume the dishes to distribute homogeneously in a square. 
This results in a homogeneous u-v sampling and alleviates 
the work involved in map making. This allows us to fol- 
low the usual cosmic shear measurement strategy and work 
in the real space. Real survey configurations often result 
in inhomogeneous u-v sampling a nd induced sy stematics in 
map making inevitably. However, IChang et alj ([2004j show 
that these systematics can be corrected. Furthermore, the 
shapelet method, which directly works in t he u-v space, is 
shown to be successful (IChang et al.|[20o3 ) and should be 
applied for real data. For these reasons, we neglect these 
systematics in our estimation. 

Since the measurement noise is often large, we do 
not try tomography in frequency space for each 21cm 
absorption "halos". We average the measurements along 
the line of sight over each "halo" and define the aver- 
aged brightness temperature along the line of sight r21 = 
^{dI/dT)5Tdv/[dI/dTAu] ~ / STdu/Av. Here, I is the 
CMB intensity and {dI/dT)5T is the associated CMB in- 
tensity decrement caused by 21cm absorption. The estima- 
tion shown here is by no means optimal and S/N could be 
improved by more advanced methods. 

The two orthogonal components of the ellipticity e of a 
21cm absorption "halo" are given by the estimator 

Q11-Q22 „ 2Q12 



el 



e2 



Qll + 022' Qll + 022 

The quadrapole moments Qij are defined by 

/ qT[T2i{e)]eiejcPe 



Qij 



J qT[T21{e)](Pe 



(2) 



(3) 



where gT(r21) is a weighting function. We choose 
qT{T21) = T21 when r21 < Tth and qT = otherwise, 
and Tth is a threshold of brightness temperature that de- 
fines the boundary of each region used for shear measure- 
ment. The position 6 is measured with respect to the center 
with i,j = 1,2 the two orthogonal directions. For a round 
object, e is induced by the cosmic shear and we have 

_ 27 27 _ 



1-K-f 7V(1-k) 



27 



(4) 



Here, 7 and k are the lensing shear and convergence, re- 
spectively. With the presence of the ellipticity eN induced 
by other sources (see § 5), we have e ~ eA'^ -I- 27. In equa- 
tion we have approximated the reduced shear <; ~ 27 
since in general k ^ 1. The second order ter m K-y in q has a 
1-10% contribution to the power spectrum (jPodelson et alJ 
|2006| ). which slightly improves the detectability. However, 
given the uncertainties in eA'', we neglect this higher order 
correction. 

An implicit assumption of the estimator (eq. 3} is that 
cosmic shear does not vary across the integral area. This as- 
sumption holds closely for galaxies, which are of arc-second 
size. The 21cm absorption "halos" are of arc-minute size and 
the cosmic shear does vary across these regions. In this case, 
the estimator (eq. |4| measures the averaged cosmic shear. 
Higher order moments can be e xplored to extract the gra - 
dient of the cosmic shear (e.g.. iGoldberg fc BaconI |2005| ). 
The relatively large size of 21cm absorption "halos" could 
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improve the S/N of such measurements. In this paper, we 
only discuss the mean shear signal of each 21cm absorption 
"halo". 



5 STATISTICAL ERRORS 

Both the system temperature noise and inhomogeneities in 
neutral hydrogen densitjjf] and spin temperature induce ef- 
fective ellipticity eN to 21cm absorption regions. This eA'' 
has expectation value zero and does not correlate over large 
distance. Analogous to galaxy intrinsic ellipticity, it only 
induces shot noise to the shear measurement. Since inhomo- 
geneities and anisotropies in spin temperature, if exist, are 
mainly caused by the central source, they are uncorrelated 
with inhomogeneities in neutral hydrogen density. The r.m.s. 
fluctuation in eA'' then has three independent contributions: 



(5) 



where (e^)|noise, (e^)|5_ff, and {e^)\Ts are ellipticity fluc- 
tuations caused by the system temperature noise, the in- 
homogeneity in the neutral hydrogen distribution, and the 
anisotropic distribution of the spin temperature around the 
first galaxies, respectively. 



5.1 Ellipticity induced by the noise in system 
temperature 

The distribution of the system temperature noise is Gaus- 
sian with an r.m.s. aT given by equation ([1]). To estimate 
the fluctuation (e'^)|noise in the ellipticity induced by this 
noise, we generate a series of realizations of observations 
by adding Gaussian noise with an r.m.s. fluctuation aT to 
the otherwise spherically distributed signal and measure the 
ellipticity using equation ([2]). The result is shown in the bot- 
tom left panel of Figure (2] For aT ~ lOmK and resolution 
0.3', the induced ellipticity is ~ 10% for z = 25. To have 
(e^) jnoise^''^ < 0.1, aT < 7 mK is required. This translates 
to a total collecting area > 300 km^ for L = 40 km. The 
induced ellipticities at other redshifts (e.g., z = 20 or 30) 
are similar, due to their similar signals (Fig. [2} . 

We notice that even for aT — > 0, (e^)|noise does not van- 
ish, rather it reaches a plateau instead. This is caused by the 
limited spatial resolution. With a finite resolution, one can 
not exactly identify the physical center of each regions, and 
any shift from the physical center would result in an effec- 
tive ellipticity. As shown in the bottom left panel of Figure[2l 
with fixed aT, one can indeed reduce {e'^)| noise significantly 
by improving the angular resolution 9p. However, in reality, 
the requirement to improve resolution is often in conflict 
with the requirement to reduce the noise. To improve the 
angular resolution, longer base line is required {9p oc L~^). 
If the total collecting area is flxed, /cover would decrease 
(/cover oc L~^). This would cause the noise aT (per reso- 
lution pixel) to increase {aT oc /"^cover oc oc 9p~^) and 
(e^) I noise to increase as well (Fig. [2)|. Therefore, improv- 
ing angular resolution does not necessarily reduce errors in 



^ Peculiar velocities are highly coherent across the 21cm absorp- 
tion region, so it is unlikely to induce a non-negligible ellipticity. 



shear measurement, unless the total collecting area is in- 
creased accordingly. Optimal surveys may balance between 
the angular resolution and aT to reduce (e^)|noise^^^ to a 
level well below the intrinsic ellipticity induced by the spin 
temperature anisotropy (fj5]3}. 



5.2 Ellipticity induced by density inhomogeneities 

The fluctuation {e^)\5H in the ellipticity induced by inho- 
mogeneities in neutral hydrogen density can be estimated 

as 



qTe'^<fk 



j qT{T21)qT{T2l')w{e - &) 



'.[{61^ - e2'^)(ei''^ - 62''^) + 4919291' e'2]d'^ed^ 



Here, w{9) is the angular correlation function of the neutral 
hydrogen over-density. In this expression, we have neglected 
fluctuations in the denominator Q11 + Q22 (eq.[2]) since they 
only cause second order effect, due to non- vanishing Qll -|- 
Q22. For a flat w{9), (e^)\5H vanishes due to the geometry 
terms in equation ((Gj. In reality, w{9) varies slowly. So w{9) 
does not contribute much to the integral in equation (|6]), 
although itself can reach ~ 10~^ at sub arc-minute scale. 
For a 21cm absorption "halo" of a galaxy in a lO*Af0 halo 

at z = 25, {€^}\5H < lO"* < ^j^) for Tth 50mK. 

Therefore, this source of error seems to be negligible. 

First galaxies are more likely to reside in high density 
regions and one may think that this strong formation bias 
invalidates the above conclusion, which is based on the en- 
semble average of the correlation function. As the shear mea- 
surement only depends on the variation in the local density 
field, to account for the environmental effect, one should 
replace ■w{9) in equation © with the conditional correla- 
tion function w{9\{5H)), where (SH) is the overdensity av- 
eraged over the 21cm absorption regions for shear measure- 
ment. Being several Mpc/h in size, these regions are much 
larger than the nonlinear scale or filament size at z ~ 25, 
which implies that (SH) <^ 1. Hence we do not expect that 
w{9\{SH)) > w{9). On the other hand, even if w{9\{SH)) 
could be 10 times larger than w{9), the induced {e^)5H^^^ 
would be still less than 3%. Based on these arguments, we 
conclude that the ellipticity induced by inhomogeneities in 
neutral hydrogen density, if not negligible, is at most sub- 
dominant to other error sources. 



5.3 Ellipticity induced by spin temperature 
inhomogeneities 

The fluctuation {e^)\Ts in the ellipticity induced by 
anisotropies in the spin temperature distribution is most 
uncertain. The spin temperature is a weighted average of 
the CMB temperature and the gas kinetic temperature. For 
the temperature range in our study, the CMB temperature 
is perfectly homogeneous. The coupling coefficients of the 
spin temperature to the gas temperature are determined by 
the spin-ch anging collisional r a te and the Lyman-a scat- 
tering rate l|Wouthuvsenlll952l : iFieldl Il959l) . On scales we 
consider, the collisional rate should be isotropic because of 
isotropic distributions of the gas temperature and density. 
However, the spatial distribution of the Lyman-a scattering 
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rate highly depends on the origin and propagation of Lyman- 
Q photons, which could introduce considerable anisotropies 
in the spi n tem perature distribution. 

.Cenl([ 20061 ) considers Lyman-Q photons originated from 
the galaxy at the center of each halo. Photons blueward 
Lyman-a from the continuum of Population III stars in the 
galaxy redshift to the Lyman-a line-center frequency be- 
cause of the Hubble expansion. At the radius where a blue 
photon redshifts to Lyman-a frequency, the photon is con- 
tinuously scattered by neutral hydrogen atoms in the IGM 
and has little spatial diffusion. The photon would not fly 
freely until it encounters an atom that has the right velocity 
to make a large jump in its frequency. The escape of pho- 
tons slightly blueward the Lyman-a line-center frequency 
from the central galaxy is affected by the neutral hydrogen 
distribution near the center. If the neutral hydrogen column 
density in the innermost Ikpc of the halo is above lO^^cm'^ 
(which is highly possible) and the column density distribu- 
tion is not isotropic (e.g., caused by the geometry of the 
star forming region or the galactic wind), it would lead to 
anisotropic distributions of the Lyman-a scattering rate and 
the spin temperature at large radii. 

However, before a photon reaches the radius where it 
redshifts to Lyman-a and becomes strongly scattered, there 
is a probability for it to be scattered and thus change its 
propagation direction. The place where the photon encoun- 
ters a tremendous number of scatterings with little spatial 
diffusion is then no longer along the initial direction that 
the photon escapes from the central galaxy. Therefore, a 
small number of scatterings encountered by photons before 
they reach the region where strong scatterings happen tend 
to make the distribution of the overall scattering rate more 
isotropic than the initial angular distribution of the photons 
from the central galaxy. 

To quantify such an effect of isotropization, we simulate 
the Lyman-a scat tering process around the galax y using a 
Monte Carlo code (|Zheng fc Miralda-Escudil2002l '). Photons 
near Lyman-a frequency are launched from the center along 
a cone with a given open angle to represent the anisotropic 
initial distribution. Every scattering of a photon is followed 
until the photon escapes to infinity. The scattering rate at 
each position in the IGM around the galaxy is then obtained 
with its normalization set by the lumi nosity of th e central 
galaxy around Lyman-a frequency (see 'Cen"200d). We cal- 
culate the resultant spin temperature distribution and find 
that with a threshold Tth = — 50mK of the 21cm absorp- 
tion region, the ellipticity fluctuation {e^)\Ts^^'^ introduced 
by the initial anisotropic Lyman-a emission are 0.21, 0.20, 
0.16, 0.14, 0.08, and 0.03 for open angles of 5°, 15°, 30°, 
45°, 60°, and 75°, respectively. We may expect that galaxy 
formation at such high redshifts is quite irregular and thus 
the formation of a well defined disk is unlikely. If this is the 
case, it may be unlikely that the effective open angles are 
much smaller than ~30°. But even in such extreme cases, 
a small number of initial scatterings are quite efficient in 
isotropizing the final scattering rate distribution by mak- 
ing most resonant scatterings happen at a position along a 
direction deviating from the initial one. 

The situation may be even more optimistic than the 
above estimates. Besides photons originally emitted slightly 
blueward the Lyman-a frequency that redshift to the 
Lyman-a frequency, there are other sources of Lyman-a 



photons that can contribute to the pumping rate of the 
21cm line. Photons originally emitted between Lyman-7 
and Lyman limit can redshift into one of the higher Ly- 
man series resonance, and, after a few scatterings, cascade 
to locally produce Lyman-a photons near the line center 
(|Chuzhov et al.ll2006l : [Chuzhov fc Shapiro! l2006l ). Owing to 
the lower cross-section of the photons between Lyman-7 s-nd 
Lyman limit, they are expected to escape from the cen- 
tral galaxy more isotropically than those slightly blueward 
Lyman-a. The small number of scatterings before they cas- 
cade to produce Lyman-a photons would further isotropize 
the distribution. An other source of Lyman-a photon s is from 
soft X-ray photons. IChen fc Miralda-Escudef (|2006l ) discuss 
21cm absorption halos around first stars and argue that soft 
X-ray photons from Population HI stars play a significant 
role in (coUisionally) generating Lyman-a photons. These 
Lyman-a photons produced locally in the inner IGM region 
surrounding the galaxy are likely to be largely isotropic be- 
cause of the largely isotropic escape of soft X-ray photons 
from the galaxy. 

Based on the above investigation, we conclude that the 
ellipticity induced by the anisotropies in the spin tempera- 
ture is likely at the level of ^ 0.1 but may be much smaller. 
Depending on the nature of these objects and observation 
configurations, either the noise induced ellipticity or the spin 
temperature induced ellipticity can dominate the error bud- 
get of shear measurement. The overall induced ellipticity is 
thus expected to be ^ 10% but could be smaller. 



6 LENSING MEASUREMENTS 

The lensing signal on arc-minute scales at 2; = 25 is 2k ~ 
10% (lower right panel. Fig. [5]). From the previous section, 
the expected noise is likely < 10%. We then expect that 
the S/N of the shear measurement for each 21cm absorp- 
tion "halo" is ~ 1. This S/N is impressive, comparing to 
S/N~ 1/30 of the conventional shear measurement of opti- 
cal galaxies a,t z — 1. 

Combining cosmic shear measurements of all "halos", 
one can measure the lensing power spectrum CI, where / is 
the multipole. The statistical error in CI, assuming Gaus- 
sianity, is 

Here, the first term is the cosmic variance and the second 
term is the shot noise, /sky is the sky coverage and Al is 
the size of the I bin, Ng is the number of 21cm absorption 
"halos". Since we measure the shear averaged over the size 
of each 21cm absorption "halo", one should replace CI with 
the smoothed CI = CIW^0{1), with WO being the Fourier 
transform of the window function describing each 21cm ab- 
sorption region used for lensing measurement. 

For a lower mass cut 7 x 1O^M0, there are about 10® 
first galaxies and thus 10® 21cm absorption "halos" across 
the sky at 24 < 2 < 26. The shot noise power spectrum is 
proportional to the square of = ae/Ng^^^. With ae ~ 0.1 
and Ng ~ 10®, the typical value of 77 is 10"'', for which 
precision measurement can be done up to Z ~ 5000 (several 
arc- minute scale; Fig|3]). 

The lensing reconstruction from 21cm "halos" can be 
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Figure 3. Error forecast for lensing power spectrum from shear 
measurement of z ~25 21cm absorption "halos" . The solid curves 
are the lensing power spectra (without smoothing) at 2 = 25 for 
cr8 = 0.77 and (t8 = 0.9, respectively. Dashed lines are statisti- 
cal error caused by shot noises as a function of 77 = ae/Ng^^'^ , 
which is ~ 10""* for typical values of tre ~ 0.1 and Ng ~ lO'^ with 
AZ = 0.1/ adopted. The quadrupole moment shear estimator mea- 
sures the mean shear averaged over an aperture and allows the 
reconstruction of lensing at Z < 5000 (dot lines, with smoothing 
area 4 arcmin^). Small scale lensing power I > 5000 can not be 
measured by the quadrupole moments shown here. However, it 
can be recovered by higher order moments such as the octopole 
moments. Eventually, one can recover all lensing information up 
to the limit of shot noise. 



further improved. Figure [3] shows that there is significant 
power at small scales (Z > 5000). Since the quadrupole mo- 
ment shear estimator only measures the averaged shear, the 
smoothing effect caused by the arcminute size of the smooth- 
ing area does not allow the extraction of small scale lensing 
information. By measuring higher order moments, however, 
smaller scale lensing information can be recovered. In the lit- 
erature, the octopole moments have been proposed to mea- 
sure the gradients of shear, or eq uivalently, the local shear 
(e.g.. lGoldberg fc NataraianlbOOa ). The relative large size of 
21cm absorption "halos" will make this measurement much 
easier than that of 2 ~ 1 galaxies. Eventually, by measur- 
ing all necessary moments, one can recover all lensing in- 
formation up to the limit of shot noise. This will allow the 
measurement of the lensing signal to I ~ 10* (Fig. [3)l . 

In the above estimation, we have assumed that intrin- 
sic ellipticities do not correlate over relevant scales and thus 
only contribute to shot noise. This is obviously true for that 
induced by local processes discussed in previous sections, 
such cis local density inhomogeneities and spin temperature 
anisotropy caused by the central sources. Could a fluctu- 
ating large scale field coupled to the brightness tempera- 
ture, such as the fluctuating soft X-ray background (e.g. 



iPritchard &l Furlanettoir2007l '). invalidate the above assump- 
tion? The answer is probably no. The modes capable of in- 
ducing ellipticities must be smaller than the ~ Mpc size of 
21cm absorption "halos". On the other hand, the modes ca- 
pable of inducing correlations in intrinsic ellipticities are of 
size ~ 2 X 10''// Mpc. Thus a fluctuating background 
can only induce correlations in intrinsic ellipticities for those 
modes of / > 10'', which are of little relevance to this paper. 

In the directions along the line of sight and perpen- 
dicular to the line of sight, 21cm absorption halos may ap- 
pear different, due to the evolution of these halos and the 
light cone effect, the redshift distortion, the beam and fore- 
ground (which are frequency and thus redshift dependent), 
etc. However, these physics do not induce asymmetry in the 
2D plane perpendicular to the line of sight. Since we only 
use the shape distortion in the same 2D plane to measure 
cosmic shear, these physics do not induce systematics in the 
shear measurement. 

Thus, with reasonable estimation of the 21cm ab- 
sorption "halo" signal at high redshifts, rather conserva- 
tive estimates for the intrinsic ellipticities caused by initial 
anisotropic Lyman-a emission from galaxies and reasonable 
estimates on the level of foreground induced ellipticity, the 
weak lensing application of these "halos" are very promis- 
ing. The combined lensing S/N of lO'' 21cm "halos" at 
z = 24 — 26 is comparable to those with traditional cos- 
mic shear measurements, with cosmic magnification of 21cm 
emitting galaxies, and with CMB lensing and 21cm back- 
ground lensing. 



7 SUMMARY AND DISCUSSION 

In this paper, we investigate the potential of using 21cm ab- 
sorption "halos" of first galaxies at z ~ 25 as background 
sources for gravitational weak lensing reconstruction. We 
show that the accuracy obtained using this method can be 
comparable, and may be superior, to some of the conven- 
tional methods. It is potentially very rewarding scientifically 
to detect and survey the first galaxies using future 21cm ra- 
dio experiments for this and other applications. 

There are several major uncertainties that may affect 
the results in this paper. The first one is the star forma- 
tion efficiency c* in large halos at z ~ 25. We have adopted 
c* = 0.2 in the calculation. A smaller c* would reduce the 
signal and size of the 21cm "halos" and thus make the ob- 
servations and shear measurement more difficult. We have 
checked that adopting c* = 0.1 doubles the r.m.s. of the 
ellipticity induced by the system noise for the same noise 
level (jT = 7 mK. Secondly, there are uncertainties in the 
hard X-ray heating, which depends on the energy extraction 
efficiency from black holes and the fraction of the released 
energy in the form of hard X-rays (see Cen 200^ for more 
details). If the real values are lower than what are adopted 
in our calculation, the X-ray heating would be less impor- 
tant and the 21cm absorption "halos" may form at redshifts 
lower than z = 25. As a consequence, we would have much 
more observable 21cm absorption "halos" and the statisti- 
cal errors of the lensing measurement would be significantly 
improved. If the X-ray heating is more efficient than we as- 
sume, the detection of 21cm absorption "halos" would be 
more demanding and it becomes more difficult to use these 
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"halos" for the lensing application. To reduce the above two 
types of uncertainties, we have to advance our understand- 
ing of the formation of the first objects. On the other hand, 
the observation (or even null observation) of the 21cm ab- 
sorption "halos" would lead to valuable constraints on star 
formation and X-ray heating at z ~ 20 — 30. Lastly, we have 
neglected any errors in the shear measurement induced by 
map m aking. Although in principle these errors can be cor- 
rected ijChang et al.ll2004l ). it is not clear how well the cor- 
rection is for unprecedented radio arrays required for 21cm 
absorption "halo" observations. The discussion of such resid- 
ual errors is certainly beyond the scope of this paper. 

While we focus here on the weak lensing application of 
the 21cm absorption "halos", because of their arc-minute 
size and high redshifts, these "halos" also have interest- 
ing st rong lensing signatures, which is discussed in lLi et al] 
ll2007t) . 
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